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Abstract 

We study the perturbative Odderon in the quasidifFractive process 'yM'yM — + r\ c °r\ c - 
At high energies this process is dominated by Odderon exchange and can be viewed 
as the theoretically cleanest test of the perturbative Odderon. We calculate the 
differential and total cross section, as well as the dependence on the energy and on 
the photon virtualities taking into account the effects of resummation of logarithms 
of the energy. The results are compared with those obtained with a simple exchange 
of three noninteracting gluons. We present the expected cross section for this 
process at a future Linear Collider and discuss implications for other processes 
involving the perturbative Odderon. 
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1 Introduction 



The Odderon is the partner of the Pomeron carrying negative charge parity quantum 
number. In high energy scattering processes it gives the leading contribution to pro- 
cesses in which negative C-parity is exchanged in the t-channel. 

After the concept of the Odderon had been proposed in PQ, it was for a long time 
almost exclusively discussed in the context of elastic or inclusive processes. These 
have the disadvantage that the Odderon gives only one of many contributions to the 
scattering amplitude and a clean identification of the Odderon is rather difficult. The 
only experimental evidence of the Odderon so far has been found as a small difference 
between the differential cross sections of elastic proton-proton and antiproton-proton 
scattering at the CERN ISR 2 a . Due to the low statistics of the data and the difficulty 
of extracting the Odderon contribution, however, it is not possible to interpret this as an 
unambiguous signal of the Odderon. For a more detailed review of the phenomenological 
and theoretical status of the Odderon we refer the reader to [3] . 

Recently an important change of direction in the search for the Odderon has taken 
place. Now the search concentrates on processes in which it basically gives the only 
contribution to the cross section. The cross section for such processes is in general 
smaller than for elastic or inclusive processes, but here already the observation of 
the process as such would establish the existence of the Odderon. Examples of such 
exclusive processes are the double-diffractive production of vector mesons in proton- 
(anti)proton scattering [3] or the diffractive production of pseudoscalar or tensor mesons 
in electron-proton scattering 0-^3- In ah of these processes Odderon exchange gives 
the main contribution to the cross section at high energies. Other possible contributions 
can only arise due to photon or reggeon exchange, but both of these contributions 
are under good theoretical control. Another interesting possibility is to study the 
interference between Pomeron and Odderon exchange. This is possible in the diffractive 
production of final states that can be produced both in a C = +1 and in a C = — 1 
state like for example a pair of charged pions. The interference term between the two 
corresponding production mechanisms can be isolated in suitable asymmetries like for 
example the charge or spin asymmetry. Asymmetries of this kind have been studied in 
|16|-|20|. Also here already the experimental observation of an asymmetry could firmly 
establish the existence of the Odderon. 

The first experimental search for one of these exclusive processes was performed for 
the case of diffractive pion photoproduction in ep scattering at HERA in This 
process is the one for which the largest cross section is expected, but its theoretical 
description obviously has to rely on nonperturbative techniques. Such a calculation 
was performed in making use of the stochastic vacuum model |221 1231 124j in the 
framework of the functional approach to high energy scattering developed in In 
|21j the experimental results have been compared to the expectations based on that 
calculation, and no signal of Odderon exchange has been found. The failure of the 
theoretical prediction for this process is currently not understood. 

In order to avoid the large theoretical uncertainties of nonperturbative calculations 
in diffractive pion production one can consider the diffractive production of heavy 
pseudoscalar or tensor mesons. In that case the large mass of the meson provides 
a hard scale, and one can hope that perturbation theory is applicable even for real 
photons. Here in particular the production of r\ c mesons has been considered, see 
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|H1 EH . The expected cross section for that process was in the range of several tens 
of picobarns. In a study of elastic pp scattering it has subsequently been found that 
the choice of parameters in the Odderon-proton coupling in those calculations was very 
optimistic and a realistic estimate of the cross section should be even smaller by 
at least an order of magnitude. Due to the small cross section, the process is not of 
immediate phenomenological interest, but it has turned out to be quite interesting from 
a theoretical perspective. 

That interest is related to the occurrence of large logarithms of the energy in the 
perturbative series. In the simplest possible perturbative picture the exchange of an 
Odderon is described by the exchange of three noninteracting gluons in a symmetric 
color state. In higher orders in perturbation theory large logarithms of the energy can 
compensate the smallness of the strong coupling constant, a s logs ~ 1, and one needs 
to resum these logarithms. For the case of the Odderon this leads to the generalized 
leading logarithmic approximation (GLLA) which is encoded in the Bartels-Kwieciriski- 
Praszalowicz (BKP) equation. Recently two different solutions of this equation have 
been found explicitly in t29, and in |.S()j . The former solution does not couple to the 
jr] c impact factor in leading order and is hence not relevant for the production of r\ c 
mesons. The latter solution, the so-called Bartels-Lipatov-Vacca (BLV) solution, on 
the other hand does couple to that impact factor. Its intercept exactly equals one, and 
it hence leads to a cross section which is constant with the energy up to logarithmic 
corrections. The BLV solution was recently also found in the dipole picture of high 
energy scattering [3T] . 

Although the intercept of the BLV solution is equal to the intercept of the simple 
three-gluon exchange model for the Odderon it has quite different properties. So far 
the phenomenological consequences of using the BLV solution in the scattering ampli- 
tude have been considered only in the diffractive production of r\ c mesons in [T31 Hi] . 
Interestingly, in it was found that for real photons the resulting cross section is by 
about a factor of five larger than the one obtained in |H] by using a simple three-gluon 
exchange for describing the Odderon. It is a very interesting question whether that en- 
hancement is a general property of the BLV solution or whether and how strongly it 
depends on the couplings of the Odderon to the proton and to the ~jr\ c impact factor. 
It is one of the aims of the present paper to address this question. 

All of the processes mentioned above involve the coupling of the Odderon to the 
proton. This coupling is known to be rather sensitive to the internal structure of the 
proton and it is therefore possible that due to nonperturbative effects this coupling 
is small. In that case it could be quite difficult to find the Odderon in these processes. 
It is therefore interesting to study also processes which do not involve the uncertainties 
of the Odderon-proton coupling. From a theoretical point of view the quasidiffractive 
process 7W7W — > r\ c r) c is the cleanest possible probe of the Odderon. Due to the 
large mass of the charm quark the coupling of the Odderon to the 77/ c impact factor 
can be calculated perturbatively even for small photon virtualities. Again, already the 
observation of this process at high energies would firmly establish the existence of the 
Odderon. More generally, one can study the quasidiffractive processes 7M7M — » MM 
and 7W7M — » MX with M being a heavy pseudoscalar or tensor meson. Also these 
can occur at high energies only due to Odderon exchange. Such processes have first 
been studied in [321 EH] and more recently for the case of r] c meson production in 
34 . In these studies the Odderon has been modeled as a simple exchange of three 



2 



noninter acting gluons. 

In the present paper we study in detail the process — > r\ c r\ c at high energies. 
In particular, we take into account the effects of resumming large logarithms of the 
energy in perturbation theory by using the BLV solution of the BKP equation. That 
allows us to perform a detailed study of the properties of the BLV Odderon solution 
in a completely perturbative process, that is in a clean theoretical setting which does 
not involve model assumptions about the impact factors. The properties of the BLV 
Odderon can easily be compared to those of an Odderon modeled by the exchange 
of three noninteracting gluons. Furthermore, by comparing the behavior of the BLV 
solution in this process and in the process "fp — ► r\ c p we can draw some conclusions 
about the possible origin of the enhancement obtained in the latter process for the BLV 
solution as compared to simple three-gluon exchange. Another important motivation 
for the present study is to estimate the chances of finding the Odderon in the process 
j(*)ry(*) _> r] c r] c in e + e~ scattering at a future Linear Collider. 

In section we provide the cross section formulae for the process — > r\ c r\ c . 
In particular we discuss the BLV solution and the jr] c impact factor. In section |31 
we study the resulting cross section and its dependence on the different parameters. 
After discussing some technical details of the calculation, we start with the case of 
real photons and calculate the differential and total cross sections in section l3~21 The 
applicability of the saddle point approximation for the BLV Odderon solution in this 
process is considered in 13.31 We investigate the energy dependence of the cross section 
in section 13.41 In section 13.51 we address the possibility to observe this process at a 
future Linear Collider. The case of virtual photons is studied in section 13.61 Finally, 
we discuss our results in the light of results obtained for the BLV solution in the process 
jp — > r\ c p in section 13.71 Our main results are summarized in section 0J 

2 The scattering amplitude 
2.1 High energy factorization 

We consider the process 'yM'yW — > rj c rj c at high energy and relatively small momentum 
transfer, that is s 3> \t\ in terms of Mandelstam variables. The photons in the initial 
state can both be real or virtual. The large mass of the charm quark provides a 
justification for treating the process in perturbation theory. 

At high energies the process is dominated by Odderon exchange. Diagrams involving 
quark exchange in the i-channel are suppressed by powers of the energy and can be 
neglected at the energies which we will consider below. Due to high energy factorization 
the scattering amplitude for Odderon exchange can be written in the form illustrated 
in figure ^ The amplitude is a convolution of the Odderon Green function G with two 
impact factors $ coupling the Odderon to the external particles. Symbolically, 

.4~<* U |G|* { ), (1) 

and the convolution includes the integration over the unconstrained transverse momenta 
of the gluons as will be described further below. The subscripts u and I of the impact 
factors stand for the upper and lower impact factor, respectively. For the Odderon 
Green function G one can insert either the BLV Odderon solution or the propaga- 
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Figure 1: Factorized form of the scattering amplitude for — > r\ c r\ c 



tion of three noninteracting gluons as the simplest possible model for the perturbative 
Odderon. 

In order to make the paper self-contained we collect in the following sections the 
known results for the impact factors and for the Odderon Green function, and we bring 
them into a form which can then be used to compute the above amplitude and the 
resulting cross section. 



2.2 The 7 — > r] c impact factor 

We first consider the impact factor that describes the transition of a photon (real or 
virtual) into an r\ c meson and three t-channel gluons in a color singlet state as it has 
been calculated in [7]. There it was found that the impact factor has only transverse 
components. The resulting expression reads 

*i = b j( v (k! + k 2 - k 3 )^ q' ^ i = 1 2 (2) 

V^) Q2 + 4m ' + (kl + k2 " k3)2 Q 2 + ^ + c ?)' '' 

where the sum runs over all cyclic permutations and the index i corresponds to the two 
possible transverse polarizations of the incident photon. Furthermore, q = k x + k2 + k3 
is the total momentum transfer, and we will have t = — q 2 . Q 2 is the virtuality of the 
photon, m c = 1.4 GeV the mass of the charm quark, e\ the totally antisymmetric tensor 
in two dimensions, and we have 



b=-^J^Tm Vc . (3) 
q c N c V a 

Here a s and a are the strong and electromagnetic coupling constants, respectively. The 
charm quark carries the charge q c = 2/3, and the r/ c meson has a radiative (photon) 
width of T = 7keV and a mass of = 2.98 GeV. 

The factor d abc is the totally symmetric structure constant for the color SU(iV c ) 
group. The scattering amplitude contains the contraction 

N 2 - A 40 

iabc jabc c _ ma ^_ (A\ 

d d ~ N c ~ 3 ' 1 ' 

where the last equality holds for A^ c = 3. 
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2.3 The BLV solution 

The Bartels-Lipatov-Vacca (BLV) Odderon solution &("> n ) found in [3U] is constructed 
from the known eigenfunctions E^ u ' n ' ) of the BFKL equation |351 I36j . These eigen- 
functions are labeled by a discrete quantum number n E Z, called the conformal spin, 
and a continuous quantum number v € R. The functions E^ v,n ^ were found in impact 
parameter space in |37| , and they can be obtained in transverse momentum space via a 
Fourier transformation. We will use the same symbol E^ u ' n ^ for both representations of 
the BFKL eigenfunctions. In the BLV solution the eigenfunctions of the BKP integral 
operator are constructed as 

^(kx.ka.ka) = c(u,n) (k l^ 2)2 ^' ra) (ki + k 2 ,k 3 ), (5) 

K 1 K 2 



(123) 



where the sum runs over cyclic permutations, and n needs to be an odd integer. We 
use the same normalization convention as in [1.3 j, so that the Odderon states have the 
same norm as the Pomeron eigenfunctions of which they are constructed. This leads 
to 

The functions (J5J) are eigenfunctions of the BKP integral operator with eigenvalues 

(7) 



, s N r a s 

xiy, n ) 



7T 



, 1 + n \ , / 1 + n 
2^(1) - ip — +iv) -VI — - iv 



where ip is the logarithmic derivative of the Euler gamma function T. 

The Odderon Green function in spectral representation is constructed as a super- 
position of all states with odd integer numbers n and general (real) v, 

OO 9 

G= V / due yx ^ n) \ L± 4 / ^ ^^ n )(ki,k 2 ,k 3 )^*^ TO )(k / 1 ,k^,k / 0. 

l^i \ r 2 (n-l) 2 ir 2 . (n+l) 2 i V X ' 2 ' 6> V !' 2 ' 3; 

odd n ^ "I 4 Jr "l 4 J 

(8) 

Here y = log(s/so) is the rapidity and sq is a fixed energy scale. The scale sq is 
undetermined in leading logarithmic approximation, and we will discuss possible choices 
for so in section 13.41 below. The normalization in (jHJ) is chosen in such a way that in 
the limit of vanishing coupling, a s — > 0, the Green function reduces to the exchange of 
three noninteracting gluons. 

In order to calculate the momentum integral for the scattering amplitude, we need 
to know the BFKL eigenfunctions in momentum space. In impact parameter space the 
eigenfunctions are 



\r10r20J \rior 2 o/ 



0) 



where rjj = r j — rj. In the r.h.s. we use complex coordinates for describing positions 
in the two-dimensional impact parameter space. Further we have the conformal weights 
h = (l+n)/2 + iv and h = 1 - h* = (1 - n)/2 + iv. 
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The Fourier transform of these eigenfunctions was calculated in It was found 
that the momentum space functions have the form 

E^ n )(l 1 ,h)=E^ n \h,h)+E^ n \h,l 2 ), (10) 

where E^ n ^ denotes an analytic contribution and E^? ,n ^ a part containing ^-functions. 
The analytic part reads 

E%' n \h,h) = C[X(h,h) -X(l 2 ,h)] , (11) 
where the coefficient C is 

(-i) n - - 
C = ±-±-hh(l - h)(l - h)T(l - h)T(l - h) . (12) 

The expression X can be given in terms of the hypergeometric function F(a\, a 2 ] b; z) = 
2 F 1 (a 1 ,a 2 ;b;z), 

- (I) (If 2 F i 1 ~ K 2 " ft; 2; "I) F i 1 ~ k 2 " k 2; -t 

(13) 

The two-dimensional momenta are denoted as complex numbers on the r.h.s. 

The <5-function part in (jlOj) is simpler. Denoting the total momentum transfer (in 
complex notation) by q = l\ + l 2 , it can be written as 

E^ n \l u l 2 ) = [6^(h) + (-1)^( 2 )(1 2 )] ^-^^^^V* 11 ' 1 ■ (14) 

2.4 Calculation of the scattering amplitude 

We want to calculate the scattering amplitudes 

A " = w ( * 1|GW) (15) 

for different transverse polarizations i, j of the incoming photons, where we have dis- 
tributed the constant factors as in [S]. In order to compute these expressions, we have 
to evaluate integrals over the independent transverse momenta, e.g. 

($i | = y'ci 2 ki ( i 2 k 2 ^(ki,k2;q)^ (iy ' n) (ki,k2;q). (16) 

However, this four-dimensional integral reduces to a two-dimensional one |3U| . 

(fcW".")) = _ r ^— / d 2 k<P\k,q-k)E^ n \k,q-k) = —^—(^lE^) 
c(y,n) J c(u,n) 

where the reduced impact factor is 



(17) 



(k,q-k)- Q2+4 ^ + (2k _ q)2 . (18) 
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Let us now consider the infinite sum over odd values of n in (jHJ) which needs ap- 
proximation in order to be evaluated numerically. In the full Green function (JHJ) the 
exponential factor e yx ^ v,n ' clearly is of special importance to the integrand, and an 
expansion of its argument can help us determine the dominant values of n. Expanding 
up to second order around v = yields 



, x N r .(Xs 

X{v, n) 



7T 



(19) 



For values of n other than ±1 we therefore get a constant part in the Taylor expansion 
of the argument of the exponential which grows with n. We have in fact checked 
numerically the contribution of n = 3 and find that this term is already of relative size 
~ 10~ 4 compared to the leading term. Therefore we can reduce the sum to one over 
n = ±1. 

Now we can further simplify the integral that we have to calculate numerically. The 
analytic part reduces for n = ±1 to 

E^ n=±1) (k, q - k) = ±7^(1 + " ™) [*( k ' q " k) - X(q - k, k)] . (20) 

For n = 1, lfT5)l leads to 

x(^)(k, q -k) = Qy^^y^F^i-z,^;-^ 

x F ( 1 - iv, 2 - iv\ 2; -^j^J ■ (21) 



For n = — 1, we get 



and thus obtain 



X in= -V (k, q-k)= X( n=1 ) (q -k,k), (22) 



E% n =- l) {k,q-k)= (k, q-k). (23) 

From now on we choose the coordinate axes such that the total momentum transfer 
q is in 1-direction, so that in complex notation q is purely real. As we have an integral 
over both components of the two-dimensional momentum vector k, we can replace the 
second component by its negative value, which in complex notation corresponds to 
complex conjugation. In the reduced impact factor (|18|) we have to switch the sign of 
the i = 1 component due to the e symbol. The i = 2 expression remains unchanged 
because it includes only the first component of the momentum in its numerator. For 
i 7^ j in (|15|) the coherent sum over n = ±1 thus gives two equal but opposite results. 
For parallel photon polarizations, i = j, the resulting expression in (|15j) is the same 
for n = 1 and n = — 1. Thus, we can work with the expression E^ ,n ~ l \ which we will 

denote as E^ from now on. 

Let us now turn to the 5- function part for n = ±1. For the ^-function part one also 
obtains the same result for n = +1 and n = —1, which we will denote by Eg , 

H = E ^n=±i) { } = r (2)( } _ (2)( _ } 1 J_ iv T{l-iv) q^_ 

5 5 y ' M ' L V ; KH '\ 2vr r(l + iv) q V ' 
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With our specific choice q = (q, 0) we can easily evaluate the convolution of the 5- 
function part with the impact factor. We can see from (|18j) that due to the e tensor 
the i = 1 component of the impact factor vanishes for k = or k = q, that is when 
the 5- functions in ()24|) are applied. Hence the only non- vanishing contribution to the 
convolution of the (^-function part with the impact factor is the one for i,j = 2, 

(6^\EV) = 1 1 F(1 = W) ^ (25) 
2.5 Saddle point approximation 

In diffractive production of an r/ c meson in ^yp scattering was calculated in the saddle 
point approximation (SPA). For our study this approximation will not be needed, and 
our results given below will not make use of the SPA. Nevertheless, it is interesting for 
us to study the same approximation also for our scattering process in order to discuss 
our results in the context of those of ^3]. For this purpose the v integral in (jHJ) is 
approximated by expanding the argument of the exponential and other z^-dependent 
factors in the integrand (in particular the BFKL eigenfunctions) in Taylor series. In 
lowest non- vanishing order the Lipatov characteristic function is quadratic in v, 

±1) = _2^C(3)^ 2 + 0{v A ) . (26) 

7T 

The rest of the integrand is also expanded in v. The first non- vanishing term in the 
expansion of the analytic part is of first order, but this term vanishes in the momentum 
integral because it is orthogonal to the ij c impact factor, see [T3j. The first term 
that survives this convolution is of second order. The (5-function part leads to a non- 
vanishing contribution already in zeroth order in 

45(k 1 ,k 2 ) = (^)(k 1 )-^ 2 )(k 2 ))^-i. (27) 

Therefore in ^3] only the 5-contribution is calculated. We will study in section 13.31 
below how that approximation affects our results. As we will show, the SPA turns out 
to be inadequate for our process. 



3 Numerical results 

3.1 Details of the calculation 

In the following sections we will present our numerical results for the cross section for 
sy(*)j(*) r} c r\ c and for its dependence on various parameters. In the present section 
we make several general remarks relevant for those results. 

The differential cross section is obtained from the amplitudes A lJ involving the BLV 
Odderon solution as 

1,3=1 

As we have seen in the previous section, the mixed polarization amplitudes vanish. 
Furthermore, we find numerically that the i,j = 1 contribution to the cross section is 



8 



only ~ 1% of that coming from i,j = 2. Since the numerical error of our calculation is 
also on the percent level we neglect this contribution. We hence have 



da 1 „ , o . . 

Tt ~ eS?W' < 29 > 

where A = A 22 . According to the discussion above we obtain the scattering amplitude 
A as the integral 

4-2 5 7 dfc yxM ^ 2 ^ 2 + ^ ^ (6\E M ) (6\E^)t (301 
3(2vr)4 J dVe 1/2(^ + 1) c(i/)2^^ )uW ^ )l > [6U) 

— oo 

where we have the reduced impact factor <p = 4> i=2 , and the convolution of (j) with E^ u ' 
is defined in Q17j). Further, E^ has two parts as in ((10(1 . and we have set n = +1 here 
while multiplying by 2 to take into account the contribution of n = —1 as explained in 
the previous section. 

Our results below are obtained from a numerical evaluation of the integral 1)30(1 . We 
emphasize that we compute that integral without further approximations. In particular 
we do not use the saddle point approximation for our main results. The outcome of 
using the SPA is included below only in order to discuss the applicability of that ap- 
proximation. In fact we will show that the saddle point approximation is not applicable 
to our process in the phenomenologically relevant kinematical region. 

Recall that the convolutions {(f)\E^) u ^ in the integrand of (|30|) involve only two- 
dimensional integrations, see (|17|). The reduction from a four-dimensional to a two- 
dimensional integral in these convolutions occurred due to the special structure of the 
BLV Odderon solution. The integral (|30|) involves hypergeometric functions which are 
expensive to evaluate in terms of computer time. But because of the reduction to two- 
dimensional integrations in (<j)\E^) u i it is still possible to perform the integral (130(1 
using Mathematica. 

As already pointed out in section l2~H we do not take into account contributions to 
the BLV solution with quantum number n/±l. We have in fact calculated numerically 
the n = 3 contribution for a variety of values of the parameters Q 2 , s and t and find it 
to be negligible. 

One of the most interesting questions which we want to study is how the BLV 
Odderon solution compares to the exchange of three noninteracting gluons in our pro- 
cess. The latter exchange is the simplest possible perturbative model for the Odderon. 
Technically speaking it amounts to replacing the Odderon Green function G in ((15(1 by 
three free gluon propagators, 

G 39 = 5^(k[ - k 1 )^(k^ - k 2 ) p-L^ . (31) 

K 1 K 2 K 3 

The cross section for 'yM'yM — > rj c rj c with a simple three-gluon exchange was calculated 
in [SI]. We have reproduced the results of that paper in order to compare them with 
our calculations. 

Our calculation is based on the BLV Odderon solution which results from the re- 
summation of leading logarithms. It should be pointed out that there are several 
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uncertainties which are inevitable in that approximation scheme. The first of these un- 
certainties concerns the choice of the appropriate value of the strong coupling constant. 
Strictly speaking the scale of a s is undetermined in GLLA. In all our calculations we 
use a s = a s (m 2 ) = 0.38 to allow for an easy comparison with the results from [31] and 
|13| where the same value had been chosen. It should be emphasized, however, that 
already a small change in a s implies a considerable change in the cross section. This 
is due to the simple fact that the cross section contains a factor a® already from the 
coupling of the three gluons to the impact factors. 

Another important uncertainty is the choice of the energy scale sq in y = log(s/so). 
Also this scale is, strictly speaking, undetermined in GLLA and has to be chosen as a 
typical energy scale for the process. We will discuss several possible choices in detail 
in section 13.41 It will turn out that choices which appear equally natural can lead to 
quite different results. 

Of course there is a minimal momentum transfer required for the transition from 
a real photon in the initial state to a r) c meson in the final state. That minimal 
momentum transfer t m i n can be estimated to be t m - m —m^ c /s. At the energies which 
we will consider its numerical value will be very small and will not have any quantitative 
effect visible in our figures. 

3.2 Cross section for the scattering of real photons 

For the calculations in this section we choose a center-of-mass energy = 300 GeV. 
This choice is somewhat arbitrary and our main motivation for it is that later on we 
want to compare the behavior of the BLV solution in the process 77 — ► r] c r] c to that in 
the process "fp — ► r] c p. The latter was calculated for the HERA energy in ^3], and we 
therefore use the same energy here. As the process 77 — ► f] c r] c is phenomenologically 
interesting mainly as a subprocess of electron-positron scattering for example, in an 
actual collider setup there is a continuous energy range available for the 77 scattering, 
rather than a fixed center-of-mass energy. We give an estimate on the size of the cross 
section for the process e + e~ — > r] c r] c in section 13. 5[ but in the first part we are only 
concerned with the properties of the differential cross section of the subprocess. The 
squared energy s enters in the argument of the exponential as y = log(s/so), with 
the scale factor sq on which we will comment later. For this section we use a scale 
sq = m 2 c . Our numerical results for the differential cross section for real photons 
(virtuality Q 2 = Q 2 , = 0) are shown in figure [2] together with the corresponding 
results for the noninteracting three-gluon process from 34 . 

Comparing our numerical results for the BLV solution to the three noninteracting 
gluons in the t-channel, we find a huge enhancement of the differential cross section at 
small momentum transfer. The cross section calculated with the BLV solution reaches 
a maximal value of ~ 120 fb/GeV 2 at t = t m - m ~ — 10 -3 GeV 2 , and then quickly falls to 
~ 88 fb/GeV 2 already at \t\ = 0.01 GeV 2 . The simple solution exhibits a fundamentally 
different i-dependence that can fairly well be described by an exponential decay. Its 
maximal value at t = t m i n is only ~ 5 fb/GeV 2 , so there is a maximal enhancement 
factor of about 25 that comes from the interaction of the three gluons. As the BLV 
curve has a much steeper i-dependence, the two curves intersect at \t\ « 3.5 GeV 2 . 

We have also estimated the total cross sections for the two cases. For the calculation 
involving the BLV solution we get a cross section of about ~ 59 fb, whereas the 
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2 4 6 8 10 12 14 

\t\ [GeV 2 ] 

Figure 2: Results for the differential cross section % for real photons (Q 2 = 0) and 
center-of-mass energy yfs = 300 GeV 



simple three-gluon process yields 0"t o t ~ 43 fb. In both cases, no cutoff for the integral 
was needed as the differential cross section falls off sufficiently quickly with growing \t\ 
to allow for a reasonable estimate of the contribution from large \t\. We notice that the 
minimal momentum transfer t m [ n is sufficiently small at yfs = 300 GeV so that it does 
not affect the calculation of the total cross section. 

In summary, we see that the BLV solution enhances the total cross section of 77 — > 
VcVc by a factor of about 1.5 compared to the noninteracting three-gluon calculation of 
[3~4j • The dependence of the differential cross section on the momentum transfer changes 
significantly, and the region of small momentum transfer becomes more important in 
the case of the BLV solution. 

3.3 Comparison with the saddle point approximation 

Next we want to study the reliability of the saddle point approximation for our process. 
This question is primarily of theoretical interest. However, it turns out that due to the 
hypergeometric functions in the BLV solution it is numerically extremely challenging 
to calculate processes in which the BLV Odderon solution is coupled to a proton. For 
these one has to make use of the SPA, and it is therefore interesting to study the 
reliability of that approximation. For this our process is well suited since here we can 
compare the SPA to the exact result. 

Figure shows our exact results together with the result obtained by using the 
saddle point approximation as described in section 12.51 Again we have chosen yfs = 
300 GeV. We find a maximal enhancement factor of the exact calculation over the SPA 
of order 25 at t = i m i n , but quickly shrinking with increasing \t\ to an enhancement of 
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Figure 3: Differential cross section: comparison of two different approximations of the 
Odderon wave function with the exact result 



order 2. For the total cross section the SPA calculation leads to crf^ t ~ 13 fb, hence 
underestimating the total cross section approximately by a factor 5. So we find that 
the SPA in its simplest form should not be applied to the scattering process at hand, 
at least if one is interested in small values of \t\. 

This result appears surprising at first sight, as it was found that in calculations 
involving the BFKL Pomeron the saddle point approximation typically overestimates 
the cross section, the deviation from the actual cross section usually being of the order 
of 20%, see for example |38] for the case of the total hadronic cross section in virtual 
photon collisions. We therefore find it instructive to discuss the origin of the large 
deviation and its direction in the case of the BLV solution in our process. 

As was mentioned in section 12.51 in the SPA the analytic part E^' of the BFKL 
eigenfunction is completely neglected. In figure El we also show how this affects the t- 
dependence. Going to small \t\, the curve resulting from the full calculation (including 
the analytic part) has a much steeper t-dependence than the SPA curve, leading to the 
large enhancement. We have also included in this figure the cross section obtained by 
neglecting the analytic part while calculating the v integral numerically without the 
SPA. That curve is very similar to the SPA curve in its i-dependence, but is higher 
by a factor of about 1.5. Thus we see that the crucial difference is not caused by the 
approximation of the argument of the exponential, but by the fact that the analytic 
part is neglected. 

The omission of the analytic part E^ in the SPA is due to the fact that the first 
non-vanishing term in the Taylor expansion (in v) of Ej? is of second order, whereas 
the delta function contribution E^ already has a non-vanishing zeroth order term. 
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Figure 4: Comparison of the different contributions to (<p\E^); the solid line is the 
Gauss-like function I{v) representing the momentum-independent factors of the v- 
integrand (the curve is scaled and shifted, see text) 



But as the factor y in the exponent is not very large, the analytic piece nevertheless 
contributes substantially to the v integral. Thus a numerical investigation of the in- 
dependence of the different parts of the solution is needed to obtain a clearer picture 
of the importance of the analytic part. 

In figure 0] we show the real and imaginary parts of both contributions to the 
expression (<j)\E^) u at \t\ = 0.01 GeV 2 . We have not included the constant factors in 
this calculation, so that the figure shows only the relative size of the different terms. 
We recall that, before being integrated, the expression {(j)\E^) u gets multiplied with 
the corresponding lower part and the rest of the v integrand in (|30j) . The latter is given 

by 



To better understand the significance of a particular contribution to the overall result, 
we have included in figure |1] also this momentum-independent contribution I(y) to the 
integral (solid curve). It basically gives a Gauss- like-curve with a maximum Jo at v = 0. 
It is scaled and shifted in such a way that the lower horizontal axis in the figure is the 
I(u) = level and the upper one the I{v) = Iq level (see r.h.s. of the figure). 

It now becomes clear from the figure, why the SPA cannot lead to good results 
in our calculation. The analytic part (dashed lines) gets comparable in size to the 
(5-function part (dotted) already at v ~ 0.2, where the Gauss-like factor (solid line) is 
still at about 70% of its maximal value. At larger values of v the analytic part even 
contributes dominantly to the amplitude. Compared to the approximated integrand, 
where the analytic part is neglected but the v integration is done numerically (no SPA) , 
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this gives an enhancement by a factor of 12 in the differential cross section (which can 
be understood when keeping in mind that the expression (4>\E^) gets squared when 
the result of the lower part is multiplied and again squared when the cross section is 
calculated). The analytic part falls off much faster with \t\ than the (^-function part, 
so the SPA improves with increasing momentum transfer. For values of | i | > 5GeV 2 
it reproduces the t-dependence fairly well. Nevertheless, even for larger values of \t\ 
the size of the differential cross section is significantly underestimated by the SPA. We 
find that the analytic part Ea becomes negligible at small t only for extremely large y 
above ~ 100. 

3.4 Energy and scale dependence of the cross section 

The intercept of the BLV solution equals unity, so there should be no power-like energy 
dependence of the cross section. But as the continuous quantum number u of the BLV 
solution leads to a cut in the complex angular momentum plane (instead of a simple 
pole), we expect a log~ c( -*) s dependence. This can be numerically verified by keeping 
t fixed and calculating the cross section as a function of s. Again, the results of this 
calculation can be compared to the SPA to check the significance of the latter. The 
comparison with the noninteracting three-gluon exchange process does not give any 
new insight, as there is no energy dependence in that cross section. 

The saddle point approximation gives an inverse logarithmic dependence on the 
energy, as can be easily seen when keeping in mind that s only appears in the Gaussian 
exponential exp(— yc'u 2 ) in the factor y = log(s/so) (with c' = 2N c a s /ir). Similar re- 
sults are expected for the numerical calculation if this approximation should be reason- 
able. However, as pointed out in the previous section, in the domain of the momentum 
transfer that gives the largest contribution to the total cross section (i.e. the small \t\ 
domain) the applicability of the saddle point approximation is very questionable. 

Again, we have calculated the differential cross section numerically, this time varying 
s. To the results a function of the form 

f(s) = alog- b (s/s ) (33) 

is fitted with fitting parameters a and b that depend only on i To see the change of 
the s-dependence with varying t, we have performed the calculation for \t\ = 0.01, 0.1, 1 
and 10 GeV 2 . The results for a wide range of squared center-of-mass energies s together 
with the fitted curves are shown in figure 

One can see that the expected behavior is reproduced quite well by the numerical 
results. This should however not be understood as a possibility to determine the value 
of b in (|33j) as a function of t. From the fits we get values for b ranging from 1.9 at 
\t\ = 0.01 GeV 2 to 2.4 at \t\ = lGeV 2 , but in the limit of asymptotically large s, all 
these curves eventually have to approach the saddle point approximation, because it 
does become valid at some large value of y when the Gaussian factor gets so narrow 
that only the region of small v gives a sizeable contribution to the integrand. Therefore, 
in the limit s — > oo, we know that b — > 1. 

In the case of center-of-mass energies accessible at present and planned accelerators 
(up to y ~ 10), however, the fitted curves give a reasonable description of the energy 
dependence. Again, we find the SPA to be inappropriate for the parameter range of 
small \t\ and realistic s. In order to keep the figure readable, we did not include the 
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Figure 5: s-dependence of the differential cross section for real photons for \t\ = 
0.01, 0.1, 1,10 GeV 2 (top to bottom), the points are numerical results, the lines rep- 
resent fits (see text) 



SPA curves, but again they fail completely in quantitatively reproducing the numerical 
results. This can be easily seen from the fact that the exponent of the logarithm is 
determined to be around 6 ~ 2, whereas the SPA gives an exponent 6=1. It is only at 
values of s well beyond any realistic size (y S> 100) that the saddle point approximation 
becomes acceptable already for small \t\. 

We now turn to the question of the dependence of our results on the choice of the 
scale so- For all of the previous calculations we had chosen the fixed scale sq = m 2 c in 
y = log(s/so). The scale so of the energy cannot be determined in LLA since a change 
in so is formally sub-leading in the expansion of logarithms. The numerical results, 
however, do naturally depend on the specific choice. It is therefore interesting to see 
the influence of different choices of sq on the cross section. If s, t and Q are fixed, a 
change of the scale so by a factor d clearly has the same effect as replacing s by s/d 
with fixed scale. It is then straightforward to obtain the resulting cross sections from 
the results obtained above. As long as the scale sq remains small compared to s, the 
change does not qualitatively alter the results. 

Yet, the overall magnitude of the cross section is significantly changed by a change 
in the scale factor. For example, in comparison with sq = m 2 c the choice so = 1 GeV 2 , 
which is the typical mass scale for hadronic processes, changes the result by an approx- 
imately constant factor of about 2/3. Compared to this uncertainty, the numerical and 
systematical errors in our calculations are definitely negligible. Strictly speaking, the 
numerical values of our results should only be taken as an indication for the order of 
magnitude. The appropriate choice for so and hence the absolute results could only be 
determined in a next-to-leading order calculation. The qualitative dependence on the 
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momentum transfer and on the center-of-mass energy, however, is quite stable under 
changes of the scale factor. 

As soon as we consider non-vanishing virtualities of one or both of the photons, the 
choice of the scale so becomes even more ambiguous. As the virtuality of the photon 
provides another momentum scale for the reaction, it is natural to include it into the 
scale factor. The inclusion of the virtuality in so will lead to different results for the 
dependence of the differential cross section on the virtuality of the photons, as will be 
discussed in more detail in section 13.61 

3.5 Possible realizations of r\ c meson photoproduction 

In order to relate our results to phenomenology, we want to give some estimates for the 
cross section in possible future collider setups. So far, we have been concerned with the 
process 77 — * rj c r] c with real photons at the fixed center-of-mass energy y/s = 300 GeV. 
That scale was chosen having in mind a later comparison to other works concerning 
the process 7p — > rj c p. 

A possible realization of quasidiffractive double rj c production would be the photon 
collider option at TESLA. In this section we use the definitions and numbers given 
in the TESLA design report [SHI- From an electron-positron center-of-mass beam en- 
ergy of -Ebeam = 500 GeV a beam of real photons with a maximum center-of-mass 
energy y/s = 390 GeV can be produced. However, due to the production mechanism 
of inverse Compton scattering, the resulting beam is not very narrowly peaked, but 
has a maximum at about y/s = 360 GeV with a width at half maximum of about 15 
%. The luminosity for photons with an energy in this peak region is estimated as 
L 77 = 1.1 • 10 cm -2 s . As the dependence on the momentum transfer at this energy 
does not give any new insight, we do not show a figure of the differential cross section 
(the curve looks exactly like figure |2J) • The total cross section for y/s = 360 GeV is 
Ctot ~ 55 fb. This would lead to a total number of events of the order of 10 5 in five 
years of continuous running. 

Another possibility of realizing the process 77 — ► r) c r} c is directly in electron-positron 
collisions. Here the process occurs as a subprocess in e + e~ — > e + e~r] c 7] c at high en- 
ergies. The 77 subsystem in such a collision has a continuous spectrum and we have 
to integrate over the energy fractions. In [H3] this calculation was performed in the 
equivalent photon approximation (see also |40| ^ for the noninteracting three-gluon ex- 
change process. There, it was much easier to calculate the convolution integral, as the 
simple solution does not exhibit any energy dependence. Because of the large numerical 
effort of calculating total cross sections in our approach, we use the results from the 
previous section to estimate the energy dependence of the total cross section. There 
we found fairly good fits of the form oc log 6 (s/so) with values of b around 2. Therefore, 
we approximate the total cross section by 

ototOO =^olog~ 2 f-^j , (34) 

where 00 = 5040 fb is determined from the value of otot at y/s = 300 GeV. 

The center-of-mass energy y/s in the 77 subprocess is related to the beam energy 

-^beam ^S 

s = z x z 2 El c&m , (35) 
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where z\ , z<i are the fractions of the electron and positron energies carried by the two 
photons. 

In the equivalent photon approximation the energy distribution of the photons is 
given by the flux-factor / 7 (z, Q^in' Qmax)- For untagged e ± it reads (for details see 
equations (23)-(25) in 



m = » + (!-»)• log a - !g^A i . (36) 



Z7T \ Z 772-g Z 

where we use G^ax = 30mrad. The cross section for the process e + e~ — ► e + e~r] c rj c , 
which corresponds to the collision of almost real photons, is given by: 

1 1 

J dzi j dz 2 0(s - s min ) o- to t(s)/ 7 (2 : i)/7(^2) • (37) 



Here s m j n denotes the minimal squared center-of-mass energy for the process. Again, 
Ebeam = 500 GeV is used. 

If we integrate over the complete domain of possible values for s, that is 4m 2 c < 
s < E^ eam we get a total cross section for the electron-positron scattering process of 
« 55 fb, as opposed to 3.5 fb that was obtained for the simple three-gluon exchange in 
|34j . However, with a squared center-of-mass energy s = s m i n = 4m 2 c one is clearly not 
in the high energy limit. In particular, the requirement s 3> t is not met. Therefore, 
we have performed the calculation again for a minimal squared energy of 10 s m j n and 
obtain for the e + e~ total cross section a value of <7 to t ~ 7fb. The large difference 
between these values comes about because the total photon cross section rises as one 
goes to small values of s. Compared to the value at yfs = 300 GeV the total cross 
section for yfs = s m i n is larger by a factor of 50. 

Stating it very cautiously, we estimate the total cross section of the process e + e~ — > 
e + e~ r q c r q c to be of the order of 10 fb. The planned luminosity at TESLA is 3.4 • 
10 34 cm _2 s _1 39 . The resulting order of magnitude of the number of events is similar 
as in the case of the photon collider option. 

A realistic assessment of the feasibility of a measurement of our process at a future 
Linear Collider would clearly require a more detailed study of the process including 
detector cuts and tagging efficiencies. Such a study is beyond the scope of the present 
paper. 



3.6 Virtual photon scattering 

So far we have been concerned with the case of real photon scattering. In this section we 
consider the cross section for virtual photons. For definiteness, we choose a virtuality 
Q 2 = 25 GeV 2 , again motivated by the choice in JH] • The comparison of the differential 
cross section for real and virtual photons is presented in figure H3 We would like to 
point out that now and in the following we are again considering cross sections for 
77 scattering rather than for the corresponding process in e + e~ scattering. We have 
included in the figure the cross sections for the case of two real photons, one real photon 
and one virtual photon, and two virtual photons. 

In the latter two processes there is an additional natural scale in the process that 
can be included in the scale factor sq. In addition to the choices discussed in section l3~H 
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Figure 6: Differential cross section for real and virtual (Q 2 = 25 GeV 2 ) photons for 
^fs = 300 GeV. 



this clearly gives another range of reasonable choices. For the case of two real photons 
we have used again sq = m 2 c , for the processes in which at least one virtual photon is 
involved, we have chosen so = m 2 c + Q 2 . That choice appears natural in particular in 
the light of the typical energy scales occurring in deep inelastic scattering. In the case 
of a virtual photon scattering on a real photon, one could also think of using some kind 
of average momentum scale (for example the geometric or the arithmetic mean). This 
would change the results only by a factor which is almost independent of \t\. 

In figure it can be seen that the basic dependence on the squared momentum 
transfer is qualitatively similar for all three processes, but the absolute size of the result 
is quite different. The virtuality of a photon leads to a suppression of the differential 
cross section. If both photons are virtual, the cross section is further suppressed. The 
relative enhancement of the real photon scattering over the processes including virtual 
photons decreases with growing \t\. This behavior does not come as a surprise if one 
looks at the reduced impact factor (|T5|) 

^-v« wtMti£^w - (38) 

The dominant region for the momentum integration is where the gluon momentum k 
is small. For small values of \t\ the suppression by the virtuality is therefore basically 
given by a factor ~ l/Q 2 . For a value of \t\ comparable in size to Q 2 , the suppression 
is only ~ 1/2, and if \t\ 3> Q 2 the effect of the virtuality becomes negligible. 

For the total cross section this gives a strong suppression with respect to Q 2 . For 
one virtual (Q 2 = 25 GeV 2 ) and one real photon, we get a total cross section <7tot ~ 5 fb, 
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Figure 7: Q 2 -dependence of the differential cross section for one virtual and one real 
photon at \t\ = 1 GeV 2 : comparison of BLV Odderon with two different choices for the 
scale factor sq and noninteracting three-gluon exchange 



for both being virtual o"tot ~ 1 fb (for both real, the cross section was <7tot ~ 59 fb). The 
three-gluon approximation leads to cross sections otot ~ 2 fb for one virtual photon and 
Ctot ~ 0.2 fb if both photons are virtual. Thus we see that the enhancement of the total 
cross section of the BLV Odderon over the three-gluon approximation gets amplified 
when one considers virtual photons. 

Next, we want to investigate the dependence of the differential cross section on the 
virtuality of one photon (figure EJ) when the second photon is real. For this we keep t 
fixed (\t\ = lGeV 2 ) and vary Q in one impact factor. As was already mentioned above, 
the results of the numerical calculations depend on the specific choice of so which can 
now include the virtuality Q 2 . In that case the Q 2 -dependence of the cross section will 
strongly be affected by the specific choice of sq. 

In figure [7] we plot the results for three different calculations: the approximation 
by three noninteracting gluons and the numerical BLV Odderon calculation with two 
choices for the scale factor. These are: the scale that was used in [T^], so = "i 2 c + Q 2 , 
which we will call 'full' scale factor because of the inclusion of the virtuality, and a 
'reduced' version without the virtuality of the photon, sq = m 2 c - As we have pointed 
out above, there is no way to determine in LLA which scale choice is the correct 
one. Consequently the terms 'full' and 'reduced' are not meant in the sense that the 
'full' scale factor is 'better' in any sense. We see that the reduced scale factor leads 
to a steeper slope than the full scale factor. Again, we find that the three-gluon 
approximation exhibits a yet steeper slope with respect to Q 2 . Nevertheless, all curves 
have a somewhat similar dependence on the virtuality. 
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Figure 8: Q 2 -dependence of the differential cross section for one virtual and one real 
photon, \t\ = 0.01 GeV 2 (upper points) and \t\ = 1 GeV 2 (lower points) with reduced 
scale factor sq; the points are numerical results, the curves are fits (see text) 



A remark is in order here concerning the applicability of the BKP Odderon in the 
case of a virtual photon scattering on a real one. If the virtuality of the former is 
large, there is an evolution in transverse momentum along the exchanged gluons in the 
i-channel. The BKP equation, on the other hand, takes into account only evolution in 
energy but not in transverse momentum. The situation of two largely different photon 
virtualities would hence require the inclusion of DGLAP type evolution |411 1421 143j 
and is not appropriately described by BKP evolution. This limitation of our results for 
largely different virtualities should be kept in mind when interpreting the results of the 
present section. 

Ignoring that limitation for the moment, we now want to investigate the 'asymptotic' 
Q 2 -dependence. That dependence is of theoretical rather than of phenomenological 
interest, but can be useful for gaining insight into the BLV Odderon solution. Again 
we consider the case of one virtual and one real photon. In figures |H] (for the reduced 
scale factor) and El (for the full scale factor) we plot the cross section as a function 
of Q 2 for \t\ = 0.01 and lGeV 2 . The points are our numerical results, the curves are 
certain fits on which we will comment below. 

If the Q 2 -dependence is not included in so, the result is simple. In the region m 2 <C 
Q 2 < s we find a (i-independent) scaling behavior of the cross section of ^ oc Q~ 4 
(the fitted curves in figure |HJ). This can be easily understood. The virtuality of the 
incident photon appears explicitly only in the impact factors (|38|) . as the BLV Green 
function itself does not exhibit a Q 2 -dependence. The momentum integrals {(p\E^) 
receive the dominant contribution from the region of small transverse momenta. As 
long as all momentum components are small compared to \J Q 2 , the error we make by 
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Figure 9: Q 2 -dependence of the differential cross section for one virtual and one real 
photon, \t\ = 0.01 GeV 2 (upper points) and \t\ = 1 GeV 2 (lower points) with full scale 
factor sq {Q 2 included); the points are numerical results, the curves are fits (see text) 



pulling the numerator out of the integral as a factor Q~ 2 is small. In the cross section 
this gives a contribution of |Q| -4 . We have added to the figures a fit with exactly that 
Q 2 -dependence to show the agreement with the numerical results. 

If we include a Q 2 -dependence in the scale so via the 'full' scale factor so = m 2 c +Q 2 , 
the overall dependence on Q 2 changes qualitatively. In the domain m 2 c < Q 2 < s we 
basically expect a curve that resembles the |Q| -4 result from the simple choice so = m 2 c 
multiplied by the dependence that we get from section l3~H cx log~ b ^\s / Q 2 ) . In figure 
OH we show the numerical results together with a fit representing that expectation. The 
fit fails when Q 2 approaches the boundary of the domain in which it is expected to 
apply but gives a fairly good description of the slope within that domain. The slope 
with respect to \J~Q 2 is about 3.2 for \t\ = 0.01 GeV 2 and about 3.0 for \t\ = lGeV 2 , 
indicating that the slope is not entirely independent of t. 

A remark is in order here concerning the question of the twist of the different solu- 
tions to the BKP equation, and in particular of the BLV solution. This issue has been 
addressed in |441 145] , see also [101 • The twist of the Odderon solutions is an interesting 
quantity from a theoretical point of view and crucial for a thorough understanding of 
the different solutions. It should be kept in mind, however, that at high energies (or 
small Bjorken-x) the operator product expansion becomes problematic and eventually 
breaks down, see for example [U]. Our results illustrate an additional problem of the 
phenomenological aspects of the Odderon twist. A potential measurement of the Odd- 
eron twist is likely to be obscured by the effects of different but theoretically equivalent 
choices of the energy scale sq, at least as long as one works in GLLA. 
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3.7 Comparison with photon-proton scattering 

In ^3] the process 7p — ► r] c p was studied using the BLV Odderon solution. The coupling 
of the BLV Odderon to the proton is more complicated than the jt] c impact factor, and 
here the reduction of the four-dimensional integral over two transverse momenta is 
not reduced to a two-dimensional one. One therefore has to make use of the saddle 
point approximation for a numerical calculation of the cross section. The quality of the 
saddle point approximation is difficult to assess in that process. A comparison with 
our process can be helpful in this respect, since there we can compare the SPA with 
the exact results and are hence able to trace the effects of the approximation. In our 
calculations we have used kinematical parameters similar to those of ^3], and we can 
therefore hope that our results can give us an idea of what the effects of the exact 
calculation in jp — > rj c p can be. 

Before we discuss possible implications of our results for the process 7p — * rj c p 
we want to point out some general differences between the results for that process 
found in JH] and the results found here for the 77 case. As was already found in the 
approximation of three noninteracting gluons |H] the differential cross section for the 
process 7p — > r\ c p vanishes as \t\ — ► 0. This is due to the coupling of the Odderon to 
the proton impact factor and is not expected in our scattering process. In addition, the 
calculation involving the BLV solution leads to a pronounced dip in the cross section 
at \t\ « 0.07 GeV 2 . Also this is due to the coupling of the Odderon to the proton as 
was pointed out in and such a dip is not expected in 77 — ► rj c r] c . 

In ^3] a sizable difference between the BLV cross section and the exchange of three 
noninteracting gluons was found in — > r\ c p. For real photons, the cross section 
is enhanced by a factor of 5 due to resummation, that is for the BLV solution. We find 
a much smaller enhancement in our process. In addition, when calculating our process 
in SPA the results are actually smaller than for three-gluon exchange. Therefore, it 
appears that in 7p scattering the enhancement is caused solely by the coupling of 
the Odderon to the proton. In the case of virtual photons (with Q 2 = 25 GeV 2 ) the 
7P calculation in SPA gives an enhancement of the total cross section by one order 
of magnitude. This means that there is an additional enhancement factor of 2 (as 
compared to the respective three-gluon calculation) caused by the virtuality of the 
photon. We find approximately the same additional enhancement factor in the case of 
one virtual photon in our process. 

In section 13.31 we found that for our process the saddle point approximation sig- 
nificantly underestimates the cross section. As the convolution of the Odderon wave 
function with the 777c impact factor in the case of photon-proton scattering was calcu- 
lated in saddle point approximation, too, it is plausible to expect that using the full 
numerical calculation also in that process would lead to an additional enhancement that 
is approximately the square root of the one found here, i. e. by about a factor 2 for the 
total cross section. We should point out, however, that it is by no means clear that the 
saddle point approximation also underestimates the Odderon-proton coupling. Instead, 
also a completely different behavior of the full 4-dimensional integral is conceivable. It 
would be important to study this problem in more detail and to determine at least the 
direction of the effect of the saddle point approximation in that coupling. 
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4 Conclusions and Outlook 



We have studied the quasidiffractive process 7M7M — ► rj c 7] c at high energies which 
is mediated by the exchange of an Odderon. This process is considered to be the 
theoretically cleanest probe of the Odderon in perturbative QCD since it does not 
involve the uncertainties typically associated with the coupling of the Odderon to a 
proton. We have taken into account the effects of resummation of large logarithms of 
the energy by using the BLV Odderon solution. This is the only solution of the BKP 
equation which couples in leading order to the 7^ c impact factor. 

We have investigated in detail the effect of resummation in this process by comparing 
our results to the exchange of three noninteracting gluons which is the simplest possible 
model for a perturbative Odderon. For real photons we find that resummation strongly 
enhances the differential cross section at small \t\, but leads to a faster decrease with 
increasing \t\. The total cross section is consequently only slightly enhanced due to 
resummation. The enhancement due to resummation is more pronounced when one 
considers virtual photons. We find a logarithmic decrease of the cross section with the 
energy in agreement with the intercept of the BLV solution being exactly one. We 
have discussed in detail the effects of different possible choices for the energy scale sq 
which is undetermined in leading logarithmic approximation. We have investigated 
this uncertainty and find that the cross section is rather sensitive to the choice of the 
energy scale, in particular in the case of virtual photons when the scale sq can naturally 
involve also the virtuality of the photons. 

We have estimated the expected event rates for the process 77 — ► rj c rj c at a future 
Linear Collider for e + e~ scattering as well as for a photon collider option. In both 
cases the observation of the Odderon in this process appears feasible, but more detailed 
studies accounting for detector cuts and tagging efficiencies will be required to obtain 
a conclusive assessment of this process. 

All previous phenomenological studies of the BLV Odderon solution were done for 
processes involving protons. Due to the complicated Odderon-proton coupling the nu- 
merical effort for an exact calculation of the cross section is prohibitively large in these 
cases. These processes have therefore been studied in the saddle point approximation. 
In our process a considerable simplification occurs due to the special structure of the 
777c impact factor and an exact numerical calculation becomes possible. We have used 
our exact results to test the quality of the saddle point approximation in our pro- 
cess. We find that for realistic values of the kinematical parameters the saddle point 
approximation fails and underestimates the actual cross section by about an order of 
magnitude. We have identified the origin of this large deviation and have indicated 
possible implications of our result for processes in which the BLV Odderon solution is 
coupled to a proton. 

Finally, we would like to point out that also other final states can be produced 
via Odderon exchange in quasidiffractive photon-photon scattering. An interesting 
example among them is the single-inclusive process 77 — > r] c X which can be treated in 
a similar way as the process discussed in the present paper. In that process, however, 
the situation is similar to the processes in which the BLV Odderon is coupled to a 
proton. The coupling of the BLV Odderon to the impact factor does not allow one 
to reduce the four-dimensional integral over transverse momenta to a two-dimensional 
one. An exact numerical evaluation of the integral is therefore not feasible and again one 
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has to make use of the saddle point approximation. In the light of our results, however, 
it seems likely that also here that approximation gives only a relatively poor estimate 
of the actual cross section. Despite this difficulty it would be very interesting to study 
the effects of resummation also in 77 — ► r] c X and in related processes involving tensor 
mesons in more detail since they might offer a good chance to observe the Odderon. 
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